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Abstract 

Numerical simulations are used to study how fiber supercontinuum generation seeded by picosecond 
pulses can be actively controlled through the use of input pulse modulation. By carrying out 
multiple simulations in the presence of noise, we show how tailored supercontinuum Spectra with 
increased bandwidth and improved stability can be generated using an input envelope modulation of 
appropriate frequency and depth. The results are discussed in terms of the non-linear propagation 
dynamics and pump depletion. 

PACS: 42.65.-k; 42.81.Dp 
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I. INTRODUCTION 



Following its first observation by Ranka et al. in 2000 supercontinuum (SC) generation 
in photonic crystal fiber has remained a subject of intense research. Motivated by important 
applications in precision frequency metrology, initial effort focussed on developing a clear 
physical understanding of the underlying mechanisms and noise properties of SC generation 
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seeded by femtosecond pump pulses, a regime which is now very well- understood [2|, 
Research is now shifting towards detailed studies of the dynamical properties for cases where 
spectral broadening is initiated in the so-called "long pulse" regime, using picosecond to 
nanosecond pulses, or even a continuous wave pump. In fact, under these conditions and 
pumping in the anomalous dispersion regime, the SC spectral broadening has been shown 
to be associated with a very rich range of dynamical behavior, involving spontaneous pulse 
break up due to modulation instability (MI) followed by the pro pag ation and interaction 



3a Dy tne pre 



10 



11 



12|. 



between very large numbers of ejected soliton pulses 

Recent results have extended this research even further, identifying significant links with 
other areas of nonlinear physics. Particular insight into the mechanism underlying SC fluc- 
tuations in the long pulse regime was provided by Solli et al. who used a novel real-time 
detection technique to directly quantify the statistics of picosecond SC shot-to-shot noise 
[isl . This work was significant in showing that the SC fluctuations led to the generation 
of "optical rogue waves," statistically rare extreme red-shifted Raman solitons on the long 
wavelength edge of the SC spectrum. Further numerical analysis of these fluctuations in Ref. 
[l^ showed explicitly that the rogue soliton statistics exhibit non-gaussian extreme-value 
characteristics. In fact, non-gaussian statistics due to soliton collisions in the presence of 
Raman scattering had been studied earlier in the context of multichannel fibre communica- 
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22l |. but possible links with the soliton dynamics of 



a |7|] were not explored. Related studies of a fundamental nature 
in the long pulse regime have applied results from weak turbulence theory to describe the 
initial spectral broadening of solitons," and the associated semiclassical scattering 

problem has been shown to analytically model experimentally-measured SC spectra 23 1. 

From an applications perspective, controlling SC bandwidth and stability is very im- 
portant and indeed, guidelines for broadband and low noise SC generation using sub-50 fs 
femtosecond pulses are well-known Recent studies in the femtosecond regime have also 
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shown how the generated SC spectra can be controlled by varying the input conditions in a 
more sophisticated manner, for example using dual frequency femtosecond pumping to in- 



duce cascaded four wave mixing 26|_or using femtosecond envelope modulation to influence 
the Raman soliton noise properties [27]. In the long pulse regime, modified initial conditions 
can also significantly influence the output SC characteristics and stability. For example, SoUi 
et al. numerically demonstrated a correlation between rogue soliton pulse height and a low 



amplitude localized noise burst on the leading edge of the picosecond pump pulses [13|. This 
idea was extended in Ref. [l4| where numerical simulations were again used to show that 
modified rogue soliton statistics could be observed using a ~4% intensity modulation across 
the full extent of the pulse envelope. A recent eprint by Solli et al. has reported significant 
experimental results, showing that seed modulation of SC generation at the -30 dB level can 
also introduce an effective phase transition in the SC stability 24 1. 



The dramatic effect of these modified input conditions illustrate the sensitivity of the 
initial MI propagation phase of long pulse SC generation to coherent input modulation. 
However, the studies described above have considered only the effects of weak perturbations 
over a limited parameter range, whereas it might be expected that the induced SC dynamics 
will depend very sensitively on both the frequency and amplitude of any applied modulation. 
Our objective in this paper is therefore to examine the potential parameter space of input 
modulation more extensively, focussing on how both the SC spectral intensity and stability 
are modified as a function of the modulation parameters. A major result is the discovery 
that certain modulation parameters yield stabilised spectra with Raman soliton peak power 
statistics transformed from an "L-shaped" extreme-value distribution to a near-gaussian 
distribution with significantly reduced peak power fluctuations below 5%. We also extend 
previous physical discussions of the SC generation mechanism in the long pulse regime to 
explicitly include the role of pump depletion. These results show that considerations of pump 
depletion dynamics can provide useful insight into the sensitivity of the SC broadening to 
input pulse modulation. 

II. NUMERICAL MODEL 

Our simulations use a generalized form of the well-known nonlinear Schrodinger equation 



suitable for modeling propagation of broadband unidirectional fields 
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inclusion of higher-order hnear and nonhnear terms, it can be written in the following way: 
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Here A{z,t) is the field envelope and the /3k's and 7 are the usual dispersion and non- 
linear coefficients. In our simulations we consider propagation in 25 m of fused silica 
based optical fibre with zero dispersion around 1055 nm. The dispersion coefficients at 
a pump wavelength of 1060 nm are: f32 = —0.820 ps^ km~^, (3^ = 6.87 x lO^^ps^km"^, 
/?4 = -9.29 X 10-^ps^km"\ /?5 = 2.45 x 10-^ps^km"\ /?6 = -9.79 x 10-i°ps6km"\ 

= 3.95 X 10-12 pg7i^^-i^ ^ _i_i2 X 10-i^ps^km-\ pg = 1.90 x 10-i^ps9km-\ 
PiQ = —1.51 X lO^^*^ ps^^ km~i and 7 = 0.015W-im-i. The nonlinear response R{t) = 
(1 — /R)(t) + fRhnit) with /r = 0.18 includes instantaneous and Raman contributions. The 
response is determined from the experimental fused silica Raman cross-section |29l| . but 
similar results can be obtained using analytic approximations for the Raman term jSOj • The 
envelope self-steepening timescale Tghock = 0.658 fs. Noise is included in the frequency do- 
main through a one photon per mode background, and via the term which describes 
thermally-driven spontaneous Raman scattering [s], [sil. To examine the SC stability char- 
acteristics, multiple simulations were carried out in the presence of different random noise 
seeds, allowing the spectral and temporal structure of the SC to be examined from shot-to- 
shot. In this regard, we have found that the spectral shape and coherence properties can 
be well determined using an ensemble of typically 100 realizations, although for cases when 
histogram data is calculated, obtaining sufficient data in the distribution tails requires larger 
ensembles with > 1000 realizations. 

We consider an input field of the form: A{0,T) = \^Pq exp{—t/2TQ) [1 + oq exp(— ifii)], 
which corresponds to a modulated gaussian pulse envelope with intensity modulation con- 
trast 2ao/(l + cio)- A modulated envelope of this form can be generated through the beating 
of a pulse with a (Stokes) frequency-shifted replica, and the use of such dual frequency 
fields has been previously used to generate high-repetition rate pulse trains through induced 



Ml/four wave mixing 32|, |33|, l3^. We carry out simulations over a range of values of oq and 



Q, adjusting the pump and Stokes sideband amplitude to ensure that the energy remains 
constant. With this approach, the individual peak powers of the pump and Stokes pulses 
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are Pp = Po/(l + clq) and Pg = Poal/{l + Oq) respectively. We choose an input energy of 
0.4 nJ which corresponding to Pp = Pq = 75 W at zero modulation. 



III. NUMERICAL SIMULATIONS 



Illustrative Results 



We first consider SC generation in the absence of any envelope modulation (oq = 0) seeded 
by pulses with Tq = 3 ps and Pq = 75 W. At this (anomalous dispersion) pump wavelength, 
the soliton order = (7P0 Tq/|/52|)°'^ ~ 111. Such a large value is typical of picosecond SC 
generation, and is expected to yield significant shot-to-shot variation in both the spectral and 
temporal characteristics 5|. This is illustrated explicitly in Fig. [D^a) which shows results 
from an ensemble of 1000 simulations, plotting the output SC characteristics after 25 m 
propagation. Subfigures (i) and (ii) show the output spectral and temporal characteristics 
respectively, superposing a subset of 75 realisations from individual simulations (gray traces) 
and also showing the calculated mean from the full ensemble (black line). 

Fig. 1 clearly shows significant shot-to-shot fluctuations in both the temporal and spec- 
tral intensities. Fluctuations in the spectral phase across the SC can be seen by plotting 
the spectral coherence calculated from the ensemble, and this is plotted on the right axis 
of subfigure (i). This coherence function corresponds to the fringe visibility at zero path 
difference in a Young's two source experiment ^ performed between independent SC spectra, 
and is a widely-used measure of SC stability [35]. Of significance here is that the shot-to- 
shot fluctuations essentially lead to zero spectral coherence (spectral phase stability) across 
the full bandwidth of the SC, except for a narrow wavelength range in the vicinity of the 
pump. Such instability in both the temporal and spectral properties is a well-known char- 
acteristic of SC generation using picosecond pulses, arising because spectral broadening in 
this regime is seeded from an ultrafast modulation that develops on the pulse envelope from 
noise-induced MI. As the envelope subsequently breaks up into individual soliton pulses 
with further propagation, the random nature of the initial modulation introduces significant 
shot-to-shot differences in the subsequent Raman soliton dynamics and frequency shifts. 

In this regard, the spectral plot in subfigure(i) of Fig. 1(a) illustrates the small number of 
events that are associated with large wavelength shifts and the generation of "rogue" soliton 
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FIG. 1: Simulation results for SC induced by pulses with (a) no input modulation and (b) mod- 
ulation at 5 THz and oq = 0.6. For each case, (i) and (ii) show output spectral and temporal 
characteristics respectively, superposing individual simulation results (gray traces) with the calcu- 
lated mean from the ensemble(black line), (i) also plots calculated coherence on the right axis, 
(iii) shows the peak power histogram after spectral filtering at 1180 nm, plotting normalized event 
frequency such that bar height represents the proportion of data in each bin. Note that coherence 
and statistical data is calculated from the full ensemble of 1000 realisations, but for clarity, the 
gray traces show only 75 realisations. 

These particular events can be isolated applying the technique developed in 
13| using a spectral filter to select components above a particular wavelength on the 
long-wavelength edge. Fourier transformation yields a series of ultrashort pulses of varying 
power depending on the position of the filter relative to the SC spectral structure. The 
statistical frequency distribution of the pulse peak power then readily reveals the presence 
of rare high peak-power events that have been fully captured because of their extreme long 
wavelength shifts. Subfigure (iii) shows the corresponding histogram obtained from the data 
using a filter at 1180 nm illustrating the heavy-tailed "L-shaped" probability distribution 
which is typical of extreme value processes. 

As discussed in the introduction above, the results in Ref. [l^ showed how a weak envelope 
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modulation could significantly modify the rogue soliton statistics. In this study, we have 
examined the effect of such input modulation for the case of higher modulation depths and 
have found that modulation depths above 50% in fact yield a very significant degree of 
control into the output SC properties compared to the case where the SC develops from 
noise. This is conveniently illustrated by considering a specific case, and Fig. 1(b) shows 
the output SC characteristics for an input field consisting of a pump at 1060 nm and a 
Stokes pulse of the same duration, but frequency-shifted by 5 THz (i.e: at a wavelength of 
1079 nm), and with oq = 0.6 such that Pp = 55 W and = 20 W. 

Significantly, although the peak power at the pump wavelength is reduced compared to 
the unmodulated case, Fig. 1(b) shows clearly that the effect of this envelope modulation 
leads to significantly different output characteristics for the spectral, temporal and statistical 
properties. For example, subfigure (i) shows an increased overall bandwidth with extended 
long and short wavelength edges, improved spectral coherence and significantly less variation 
in the shot-to-shot spectral structure compared to Fig. 1(a). This leads to a well-isolated 
Raman soliton peak on the long wavelength edge. The time domain characteristics in sub- 
figure (ii) also show reduced shot-to-shot variation when compared to Fig. 1(a) and clearly 
illustrate the localized temporal structure of the Raman soliton around a time coordinate of 
24 ps. However, the inset shows that although each Raman soliton has a temporal duration 
of ~ 80 fs, the residual wavelength fiuctuations result in significant temporal jitter in the 
soliton position of ~1 ps. It is this temporal jitter that leads to the near zero spectral 
coherence in subfigure (i) in the vicinity of the Raman soliton wavelength. Nonetheless, 
when compared to the spontaneously generated SC in Fig. 1(a) the input modulation yields 
a remarkable improvement in the SC stability properties. This is also seen in subfigure (iii) 
which shows the peak power histogram after long wavelength filtering at 1180 nm. Here 
we see that instead of an L-shaped distribution indicative of a small number of extreme 
red-shifted rogue events, we obtain a significantly more localized gaussian-like distribution. 
Indeed, the fiuctuation in the filtered soliton pulse peak power about the mean for these 
parameters is only 5%. 
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B. Interpretation and Discussion 



To interpret these results physically, it is useful to consider how the SC characteristics 
vary a wider range of modulation frequencies. To this end, Fig. 2 shows results for ao = 0.6 
but over a 0-20 THz span of modulation frequencies. The results show (a) the mean output 
spectra and (b) the associated wavelength-dependent degree of coherence calculated from 
the ensemble. Here we plot the results using a false color representation with wavelength as 
the horizontal axis and modulation frequency on the vertical axis. Note that these spectral 
characteristics calculated over an ensemble generalise similar single-shot noise-free results 
used in previous work [l^, and provide more realistic predictions of the spectral structure 
and stability properties that could be expected in experiments. 

From these results, we see that the output spectral structure clearly exhibits significant 
dependence on modulation frequency. Indeed, it can be seen that the particular frequency 
of the input modulation introduces a remarkable degree of control into the SC generation 
process, allowing both extension and reduction of the output bandwidth when compared 
to the unmodulated case. The dependence of the output spectral structure on modulation 
frequency is complex, but useful insight can be obtained by comparing the results seen in 
Figs. 2(a) and (b) with the frequency-dependent gain curve describing the growth of a 
Stokes wave component due to MI and Raman processes. To this end, we plot in Fig. 2(c) 
the gain curve describing the coupling between MI and Raman gain 36| , calculated for the 
un depleted pump wave of power Pp at the fiber input. Note that the curve is more complex 
than a mere superposition of the separate MI and Raman processes, but nonetheless shows 
a clear MI gain peak around 8 THz and a distinct shoulder around the 13.2 THz peak of 
Raman gain in fused silica. 

From inspection of Fig. 2(a), we can identify three broad ranges for the input mod- 
ulation frequency that present qualitatively different output characteristics. Firstly, for 
modulation frequencies less than ~ 8 THz, the output spectra extend on both the long and 
short wavelength side of the pump and exhibit soliton and dispersive wave structure. These 
characteristics illustrate the importance of soliton dynamics in this regime, arising because 
of the reshaping of the input pulse envelope into a train of high contrast soliton train during 
the initial MI phase of propagation. Significantly, we note that the values of modulation fre- 
quency detunings where we observe these dynamics are less than the frequency of maximum 
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Wavelength (nm) Wavelength {nm) Gain (m'^) 



FIG. 2: (a) Density plot of output spectral intensity after 25 m propagation for ao = 0.6 as a 
function of applied modulation frequency (vertical axis), (b) Corresponding degree of spectral 
coherence, (c) Calculated mixed Mi-Raman gain curve. 

MI gain. Although perhaps surprising, this can be understood from the fact that we are 
in a strong conversion regime, and thus frequency of maximum MI gain dynamically shifts 
towards lower values with propagation as pump depletion becomes significant. 

In fact, previous analytical studies of this process under CW excitation using a trun- 
cated sideband model to describe coherent pump-energy exchange have predicted maximum 
integrated gain with an initial modulation frequency of around half the frequency of peak 
MI gain calculated using the initial value of undepleted pump power [sT^. Although the 
dynamics of the regime considered here is more complex (involving multiple sidebands and 
Raman scattering), such an analysis is useful for the physical interpretation of these results. 
In particular, we can interpret the improved coherence properties of the generated SC seen 
in Fig. 2(b) [see also Fig. lb)] as arising from the reduced influence of spontaneous MI 
relative to coherent exchange between the generated sidebands. 

Secondly, as the modulation frequency increases into the range 8-15 THz, Fig. 2(a) 
now shows a predominantly red- shifted output SC spectrum where long wavelength soliton 
structure is no longer clearly apparent. This behaviour can be readily understood from 
the fact that the dynamically decreasing frequency shift of the MI gain curve due to pump 
depletion means that increased initial modulation frequency lies outside the region of sig- 
nificant MI gain after a very short propagation distance and thus the dynamics exhibit 
reduced coherence as reflected in Fig. 2(b). In this regime, Raman amplification plays a 
more significant role and there is a predominantly red-shifted output spectrum. A more 
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detailed analysis on a shot-to-shot basis still shows the development of red-shifted solitons 
from the Raman-amplified Stokes side-band, but these are not distinctly resolved in the 
mean spectrum. 

Finally, for frequencies exceeding ~ 15 THz, the input modulation lies outside the regime 
of significant MI and Raman gain. In this case, we see the growth of low amplitude sponta- 
neous MI sidebands around 7 THz from the pump so that this regime is essentially similar to 
the incoherent and unstable case of an unmodulated pump in Fig. 1(a). A detailed consid- 
eration of this regime would, however, need to take into account energy exchange between 
the pump and both initial and spontaneously generated sidebands with propagation. 
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FIG. 3: Energy over 1 THz bandwidth at the pump wavelength as a function of propagation 
distance. Left: false color representation for modulation frequencies over the range 1-20 THz 
up to a propagation distance of 7 m. Right: pump energy over 25 m propagation for selected 
modulation frequencies as shown. 

Fig. 3 and 4 illustrate additional aspects of the dynamics for this case of a 5 THz 
modulation and with = 0.6. Specifically, Fig. 3 highlights the importance of pump 
depletion in the initial phase of propagation, plotting the energy at the pump wavelength as 
a function of propagation distance, and for modulation frequencies over the range 1-20 THz. 
Note in calculating the energy at the pump wavelength we integrate over a 1 THz bandwidth; 
this is more than 10 times the initial bandwidth, and allows the depletion of pump energy 
through gain dynamics to be fully captured while covering a bandwidth sufficient to include 
the broadening of the pump itself due to self-phase modulation. 

The left subfigure shows results using a false color representation up to a propagation 
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distance of 7 m, whilst the right subfigure plots the pump energy over 25 m propagation 
for selected modulation frequencies as shown. Although the detailed frequency and distance 
dependence of the pump depletion is complex, it is nonetheless clear that maximum pump 
depletion occurs at an earlier stage of the propagation for modulation frequencies less than 
8 THz, correlating with coherent MI gain dynamics and the more distinct observed soliton 
features in the spectra shown in Fig. 2. 
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FIG. 4: Dynamics for the spectral broadening in Fig. 1(b) with modulation frequency of 5 THz and 
amplitude oq = 0.6. Left: line plots of temporal intensity profile at selected distances. Right: false 
colour plot of the intensity profile evolution, with a detailed view showing how the largest soliton 
component collides with and extracts energy from other soliton pulses on the pulse envelope. 



The particular modulation frequency of 5 THz considered in Fig. 1 leads to a very clear 
and distinct soliton component on the long wavelength edge of the output spectra, and Fig. 
4 considers the corresponding time-domain dynamics in more detail. This figure shows how 
the evolution of the initially modulated envelope leads to the ejection of one particular high- 
power soliton pulse (after a distance of around 7 m) that extracts energy from the other 
soliton like pulses on the envelope, before clearly separating from the residual envelope 
in the time domain due to its significantly different group velocity. As mentioned in the 



introduction, the exchange of energy between colliding so 



effect has been the subject of much previous research 
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itons in the presence of the Raman 



2l| and it is likely that this process 



plays a key role in the statistical excitation of rare rogue wave events in the case where 
supercontinuum generation arises from noise Q, [l^. The significance of the results shown 
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FIG. 5: Top: false color density plots of output spectra as a function of modulation frequency for 
four values of modulation index as indicated. Bottom: corresponding plots of the mutual coherence 
function. 

in Fig. 4 is that they suggest that the study of the propagation dynamics of an initially 
modulated input may allow this energy exchange process to be studied under controlled 
(rather than statistical) conditions. 

As a final point of discussion, we note that similar qualitative features are observed over 
a wider range of modulation depths, and Fig. 5 shows false colour density plots showing the 
modulation frequency dependence of the mean output spectra and associated coherence for 
different values of the modulation depth. Although the detailed interpretation of the results 
must be carried out on a case-by-case basis, we can clearly see that increased modulation 
depths are associated with improved coherence in the 0-8 THz range where coherent MI 
gain processes are significant. This can be readily understood from the fact that the strong 
conversion regime where coherently seeded MI initiates the SC process is only reached for 
larger values of the modulation depths, in which case the influence of spontaneous MI is 
strongly reduced and overall coherence is enhanced. 
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IV. CONCLUSIONS 



In this work we have investigated the possibihty of controUing the generation of supercon- 
tinuum in the long pulse regime through modulation of the input pulse. Our results indicate 
that depending on the value of the modulation frequency and depth the mean supercontin- 
uum spectral shape can be significantly modified. In addition, we have shown that, for large 
values of the modulation depth and modulation frequencies in a range determined by the 
initial Mi-Raman gain, SC spectra with extended bandwidth and improved stability can be 
generated. These results may find apphcation in the tailoring of broadband SC spectra for 
specific applications. 
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